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Abstract. Enhancement of nuclear pasta formation due to multi-nucleus simultaneous collision is presented 
based on time-dependent density functional calculations with periodic boundary condition. This calculation cor- 
responds to the situation with density lower than the known low-density existence limit of the nuclear pasta 
phase. In order to evaluate the contribution from three-nucleus simultaneous collisions inside the cosmic matter, 
the possibility of multi-nucleus simultaneous collisions is examined by a systematic Monte-Carlo calculation, 
and the mean free path of a nucleus is obtained. Consequently the low-density existence limit of the nuclear pasta 
phase is formed to be lower than believed up to now. 



1 Introduction 

Nuclear pasta, whose existence and role were suggested in 
association with the collapse driven supernova explosions 
10121 . is of great importance with respect to the astrophysi- 
cal synthesis of chemical elements. Indeed nuclear pasta is 
known to change the transparency rate of neutrinos II3I4I . 
which play a role of trigger in supernova explosions giving 
rise to a wide variety of chemical elements. Therefore the 
existence of nuclear pasta at low densities is expected to 
have a great impact on clarifying the historical evolution of 
chemical elements and also on revising the low-density and 
low-temperature part of the nuclear phase diagram |5 6|. In 
particular, here we are concerned with its existence in the 
outermost part of protoneutron stars (for the review of su- 
pernova mechanisms, see |7|). 

In this paper, based on three-dimensional time-dependent 
density functional calculations, a dynamical formation of 
nuclear pasta in the low-density and low-temperature sit- 
uations is demonstrated by 3-nucleus simultaneous heavy- 
ion collision. This calculation simulates the situation out- 
side of the known low-density existence limit of nuclear 
pasta phase. In order to examine the validity of this calcu- 
lation, the possibility of multi-nucleus simultaneous colli- 
sions is studied by a systematic three-dimensional Monte- 
Carlo calculation, and the mean free path of a nucleus at 
a given density and temperature is obtained. Note that for 
this kind of nuclear pasta formation, the presence or ab- 
sence of the fast charge equilibration process |8| deter- 
mines the upper boundary of nuclear pasta phase with re- 
spect to the temperature ||9l. 



star. In the following the dilute cosmic matter means the 
matter whose density is less than 10 % of the standard nu- 
clear density. 

Although the existence of nuclear pasta was suggested 
in more dense situations, its existence on the outermost 
part of protoneutron star has never been understood. Its ex- 
istence, however, is shown in this paper. 



2.2 Multiple nucleus collisions 

We are concerned with heavy-ion collisions appearing in 
the dilute cosmic matter. Here the mean free path of a nu- 
cleus is important, and multi-nucleus simultaneous colli- 
sions might play a role. Indeed, rather dense situations, 
which are not achieved on the earth, appear in the universe, 
so that the presence of multi-nucleus simultaneous colli- 
sions is expected. For reference, the density of targets in 
accelerators is roughly estimated by 



PLab = 10 'Vo 



14 



where the radii of nuclei and atoms are assumed to be 1 
m and 10"'" m, respectively. In the following the density is 
denoted by p, and the terminology of multi-nucleus simul- 
taneous collision exactly means simultaneous collisions be- 
tween more than two nuclei. 



3 Calculations of 3-nucleus simultaneous 
collisions 



2 Heavy-ion collisions in cosmic matter 

2.1 Dilute cosmic matter 

Cosmic nuclear matter, which consists of both neutrons 
and protons, is studied. Here we are concerned with the 
dilute cosmic matter whose density is less than 10 % of the 
standard nuclear density po- Nuclear matter of this kind is 
expected to exist in the outermost part of the protoneutron 



A three-nucleus simultaneous collision is shown based on 
three-dimensional time-dependent density functional cal- 
culations employing the Skyrme force parameter set SLy6 
lllOJ . The calculation is carried out in a spatial box 64 x 
32 X 64 fm^ with a spatial grid spacing of 1 .0 fm, where 
periodic boundary condition is applied. Three ^^Zn are ini- 
tially located at (8,0,-8), (-8,0,-8) and (0,0,8), respectively 
(cf. the first frame in Fig. [T]). The two nuclei initially lo- 
cated at (8,0,-8) and (-8,0,-8) are not given velocities in the 
center-of-mass frame, while the nucleus initially located 
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Fig. 1. (color online) Time evolution of a 3-nucleus simultaneous collision in a cubic space with periodic boundary condition. The box 
size is 64 fm x 32 fm x 64 fm. 



at (0,0,8) has the initial velocity along the z-axis. Accord- 
ingly, the total system consists of 90 protons and 144 neu- 
trons, and the total incident energy is set to 100 MeV in the 
center-of-mass frame. 

Figure [1] demonstrates a 3-nucleus simultaneous colli- 
sion resulting in nuclear pasta. The temperature and den- 
sity of this calculation are 0.5 MeV (see following discus- 
sion leading to Eq. ([T]i) and 0.039 po, respectively. This 
situation corresponds to the state with a density lower than 
the known low-density existence limit of nuclear pasta phase 
f5|. If we consider collisions with a higher energy (T = 
7.5 MeV), nuclear pasta cannot be formed. It corresponds 
to the energy above the upper-limit energy of fast charge 
equilibration |8|. For comparison, 2-nucleus collisions re- 
quire more dense situations to form a nuclear pasta. A cal- 
culation with exactly the same setting except for the force 
parameter set (SKI3 IfTTII ) results in the same conclusion. 
Consequently the low-density existence limit of nuclear 
pasta phase is suggested to be extended by taking into ac- 
count multi-nucleus simultaneous collision dynamics. It is 
important to note that the boundary condition is not es- 
sential. Therefore Fig. [T] should not be interpreted as the 
formation of infinitely connected matter, but simulates the 
favorable situation of very long finite spaghetti-like struc- 
ture. Pasta formation by 3-nucleus simultaneous collisions 
has been calculated for several different settings of the ini- 
tial condition |^ . 



4 Monte Carlo calculation 

Although multi-nucleus collisions are suggested to enhance 
nuclear pasta formation, it is worthless if such a situation 
is quite rare in the cosmic matter. 

Let us confirm the validity of multi-nucleus collisions 
in the cosmic matter The numbers of simultaneous colli- 
sions are calculated utilizing the Monte Carlo method. A 
cubic space with side length L fm is given. Spheres (corre- 
sponding to nucleus) with the same radius r fm are intro- 
duced in the cubic space, where those positions are chosen 



randomly (The left panel of Fig.|2]i. Accounting for the ef- 
fect from the Pauli principle (Fermionic feature) and the 
Coulomb force, a sphere overlapping with the other sphere 
is discarded from the beginning. Together with periodic 
boundary condition, thermal equilibrium of nuclear mat- 
ter is simulated. The density of a sphere is fixed to po, so 
that the density of matter is given by adjusting the values 
of L {L is fixed in actual calculations), r and the number of 
spheres. Another sphere (nucleus) with exactly the same 
feature (radius: r fm, density: po), which corresponds to 
the projectile, is coming from the outside of the cube. The 
temperature of the system is determined by the relative ve- 
locities between spheres. Here it is given by the relative 
velocity of the projectile to the cubic space. The projectile 
experiences collisions with spheres, when it goes though 
this periodic space. 

Let the typical time interval of low-energy heavy-ion 
collisions be 1000 fm/c = 33.3 x 10"^^ s. Collisions are 
defined to be simultaneous, if they occur within this time 
interval. The statistical values for the number of simulta- 
neous collisions is obtained by lO'* trial calculations. Let A 
denote the nuclear mass number. Three cases with different 
radii of the sphere (A = 50, 100, 200): 



r= 1.2xA'/^fm, 

are calculated for several densities and temperatures. In 
particular, with respect to low-energy heavy-ion collisions, 
we consider three different temperatures T = 1,5, 10 MeV. 
According to the Stefan-Boltzmann relation, the tempera- 
ture of matter in heavy-ion reactions is estimated by the 
relative kinetic energy. 



o-r = E\ 



(1) 



where T and E* mean the temperature and the incident en- 
ergy, respectively. 
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Fig. 2. (color online) In the left panel the setting of the Monte Carlo calculation is shown, where periodic boundary condition is applied. 
Spheres, whose spatial distribution is random, stand for nuclei. The passing of the projectile going through nuclear matter is simulated. In 
the right panel the result of Monte Carlo calculation is shown. The number of collisions during a typical interval of low-energy heavy-ion 
collisions is presented. 



Table 1. Mean free path [fm] of a nucleus in nuclear matter, 
which consists of nuclei with mass number A. 
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4.1 Mean free path of nucleus in the matter 



Table[T]shows the mean free path of a nucleus at given tem- 
perature and density. The mean free path tends to be larger 
for matter including heavier nuclei. It seems to be reason- 
able; for a fixed density and volume of the cubic space, 
there are more nuclei in lighter case. The mean free path 
becomes smaller for dense and high temperature situations. 
Comparing two cases in p = 0.02po and T = 1 MeV, al- 
most 35 fm difference of mean free path is noticed depend- 
ing only on the size of the nuclei. Note that the cases with 
r = 1 MeV roughly correspond to heavy-ion reactions 
at the incident energy slightly higher than the Coulomb 
barrier. This table implies that multi-nucleus simultaneous 
collisions commonly appear in the dilute cosmic matter 
satisfying p > O.Olpo. 



Table 2. Transition density [10""po] from dominance of 2- 
nucleus collisions to that of multi-nucleus collisions. The dom- 
inance of 2-nucleus collisions is defined to be true, if 2-nucleus 
collisions are the most probable among other finite-nucleus col- 
lisions (not taking into account coUisionless cases). The density 
dominated by the 2-nucleus collisions (the left side of arrow) and 
that dominated by the multi-nucleus collisions (the right side of 
arrow) are shown in each case. 
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4.2 Transition to the situation dominated by 
multi-nucleus collisions 

The right panel of Fig. |2]shows the statistics of simultane- 
ous collisions for T = 1 MeV and A = 100. It demon- 
strates the density-dependent transition from dominance 
of 2-nucleus collisions to that of multi-nucleus collisions. 
For the lowest density case with p = 0.02po, the situation 
is dominated by 2-nucIeus collisions. However, depend- 
ing on the increase of the density, 2-nucleus collisions lose 
their dominance below p - O.lOpo, and multi-nucleus col- 
lisions becomes dominant instead. It is also seen that three- 
nucleus simultaneous collisions are not expected (their ex- 
pectation values are negligibly small), if we are only inter- 
ested in collisions in the laboratory (p = pLab)- 

Table |2] summarizes the transition from the situation 
dominated by 2-nucleus collisions to that dominated by 
multi-nucleus collisions. The transition density is found to 
be of the order of lO^^po, where its temperature depen- 
dence is small. In addition simultaneous collisions between 
many nuclei are expected for the nuclear matter with 1 % 
of the standard density. Indeed, for p - po and T = 1 MeV, 
the most frequent numbers of collisions are 77, 59, 47 for 
A - 50, 100, 200, respectively. Note again that this cal- 
culation is carried out by giving a typical time interval 
(1000 fm/c), so that the transition density becomes larger if 
we assume the smaller time interval. Consequently, multi- 
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nucleus simultaneous collisions become important and are 
expected to play a role in the dilute cosmic matter with 
density higher than 0.0 Ipo, and the validity of the calcula- 
tion shown in Fig.[T]is confirmed. 



5 Summary 

Enhancement of nuclear pasta formation due to 3 -nucleus 
simultaneous collisions has been presented based on time- 
dependent density functional calculations (using SLy6 and 
SKI3 force parameter sets). Multi-nucleus simultaneous 
collisions between more than two nuclei have been con- 
firmed to occur in cosmic matter with density higher than 
O.Olpo- Therefore, nuclear pasta formation is possible in 
the outermost part of protoneutron stars, and neutrino radi- 
ation is not as easy as expected. Eventually, the low-density 
existence limit of nuclear pasta phase is suggested to be 
lower than beUeved up to now. 
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